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Abstract 

NLO QCD corrections to charged and neutral current Drell-Yan processes and their 
implementation in the computer system SANC are considered. On the partonic level 
both quark-antiquark and quark-gluon scattering channels are taken into account. Sub- 
tractions of the collinear singularities in the massive case are compared with ones in 
the MS scheme. Results of SANC on the hadronic level are presented. Comparison with 
results of the MCFM package is shown. 

1 Introduction 

Charged and neutral current Drell-Yan (DY) processes [I] on the eve of the first proton 
collisions at the LHC become very important for precision tests of the Standard Model. 
They are easily detected and will provide standard candles for detector calibration during 
the first stage of LHC running. They will be also used for extraction of partonic density 
functions (PDF) in the kinematical region which has not been accessed by earlier experiments. 
Therefore it is crucial to control the theoretical predictions for production cross sections and 
kinematic distributions of these processes. 

In the previous paper [2] we presented a part of the QCD sector of our computer system 



SANC [5] {http://sanc.jinr.ru/ and http://pcphsanc.cern.ch/ ) where the NLO QCD processes 



are treated. There we considered the implementation into SANC the calculation of the charged 
(CC) and neutral (NC) current quark-antiquark Drell-Yan processes on the partonic level 
and briefly presented some numerical results for the hadronic level. The QCD corrections to 
DY processes are known in the literature for many years, see Refs. [HI El [8]. Recently the 
corresponding NNLO corrections for differential distributions have been received P, [10] . 
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In this paper with respect to Ref. [2] we add into consideration quark-gluon and gluon- 
antiquark Drell-Yan processes on the partonic level side by side with the quark-antiquark 
ones. We implemented into SANC the calculation of QCD corrections to DY processes on the 
hadronic level. Working with massive quarks, we regularized the collinear singularities by 
masses of quarks. But on the hadronic level we have to remove these collinear singularities 
to avoid double counting, because they are already included in PDF. Therefore we compared 
analytical results obtained in our treating of collinear singularities with analogous results 
obtained in the MS scheme calculated in the n-dimensional phase space. In this way we 
extracted the subtraction terms needed to remove the collinear singularities in our massive 
quarks case. We show here also comparison with the corresponding results of the MCFM [TT] 
package. 

One-loop electroweak radiative corrections were computed for the DY processes by SANC 
in Refs. [T2j [T3 | [H] and extensively compared with results of other groups, see e.g. Refs. [T5| 

[HE]. 

The paper is organized as follows. In the second section we calculated the hard gluon 
bremsstrahlung contributions both to charged and neutral current Drell-Yan processes on a 
quark-parton level in the massless MS scheme and compared with our massive quarks results. 
We used F0RM3 . 1 [18J for analytical calculations. In the third section we calculated the quark- 
gluon and gluon-antiquark DY processes with charged and neutral currents in our massive 
quarks treating and draw parallel with calculation in the massless MS scheme. In Conclusion 
we discuss the results and give some illustrations of DY distributions at the hadronic level. 



2 Quark— ant iquark Drell-Yan processes 



2.1 Massive quarks treatment of the hard gluon contribution 

Working with massive quarks we showed in the previous paper [2] the NLO QCD corrections 
due to the hard gluon bremsstrahlung of charged current (CC) d(pi) + u{p2) — > W — > 
ve(p3)+t(p4)+g(P5) and neutral current (NC) q(px)+q(j?2) -> [A, Z] -> £{p 3 ) +I(p 4 ) + g{p<r>) 
Drell-Yan processes were obtained in the form 
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where \x is the factorization scale. The energy of the emitted gluon is 



(2) 



(3) 
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s' 

where z — — , s — — 2(pi.p 2 ), p\ and p 2 are momenta of the incoming quarks and s' = 
— (p3 + Pi) 2 is the invariant mass of the outgoing leptong^; 

P„{z) = C F ±±£ (4) 
i — z 

is the leading order (LO) quark-quark splitting function. 

2u> 

In the charged current case we have the limits of integration over variable z: z max = 1 -=., 

vs 

where uj <C \/I and z m i n = — -, mg being the mass of charged leptons. The auxiliary parameter 

s 

uj is the maximal energy of a soft gluon in the c.m.s. of the incoming partons. In the neutral 

4/77 ^ 

current case we have correspondingly the same z max and z min = -. 

s 

The cross sections in the Born approximation read 
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where s = — (pi + P2) 2 , Pi and p 2 are 4-momenta of the initial quarks; = M\ — iM w Y 
/ 4wi« 

/3(s,m 2 ) = y 1 — . Here we denoted 

W) = QjQl + 2Q ff Q< I x,(5) I v q v t + \ xz(s) I 2 (?, 2 + if 2 ) ^ + /f 2 ) , 

V a (S) = ^)-2|x^)| 2 K 2 + /f )2 )(/f) 2 , 



4 3) - 2Q, sin 2 9 W , v t = lf ] - 2Q e sin 2 9 W . (7) 



The Z/7 propagator ratio with s-dependent or constant Z- width is 

s 1 
= r~ a ■ 22 2~H~ ■ ( 8 ) 
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We see in Eqs. flTJ and fl2]) that the collinear singularities appear as quark mass singularities. 



2.2 Hard gluon contribution in MS scheme with massless quarks 

Because the collinear singularities calculated in the MS scheme with massless quarks are 
already included in PDF, we have to find which terms in our massive quarks treatment of the 
cross sections have to be subtracted to avoid the double counting. Therefore we calculate the 

1 We use the (— , +, +, +) metrics, p = {p°,p). 
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same cross sections in the MS scheme following the well known [19] manner of working with 
massless quarks and compare with our results. 

Calculating the three particle phase space element of the hard gluon emission in the n- 
dimensional phase space we used a cascade in two steps: 



ds' 



Z7T 



(9) 



For the first step of the charged current process d(pi) + u(j>2) — ► W*(Q') + g{p*) we obtained 
the same formula as in H91: 
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where we introduced in addition to z the variable y 
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dyy £ {l-yY 



(10) 
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Here 9 g is an angle between vectors p\ and p$, the angle of the emitted gluon. 

The phase space element of the second step W*{Q') — > u^{p^) + £ + (p4,) of the cascade is: 
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where 



y 



R 



1 + cos(8 R ) 



(13) 



and 9 R is an angle between the charged lepton and gluon in the rest frame of the outgoing 
leptons. 

Having in mind that the cross section of the charged current process in Born approximation 
has the form 
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we obtained a factorized expression of the hard gluon NLO correction to the charged current 
process. 
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Integration over y gives 
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One can see that the collinear divergence appears here as a pole -. To compare this expression 

£ 

with the analogous expression flTJ where the collinear divergence manifests itself in the form 
of logarithms In ( — - ] , In ( — rj I one have to take the limit e — > . Then one obtains an 



m; 



expression 
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to be compared with the corresponding expression (jTJ). We see which terms in the expression 

(PQ) correspond to the collinear divergent term — — which in MS scheme has to be subtracted 

from the hard gluon contribution to the considered process because it is already included into 
PDF. 

So, on the quark-parton level we have subtract from a^ard ©> the following expression: 
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Factorization properties and general relations between amplitudes with massive and massless 
partons can be found in Ref. [2D] . 

For the neutral current process calculating the three particle phase space element (Q of 
the hard gluon emission in the n-dimensional phase space we obtain the same result flDl) for 
the first step q(pi) + q(p2) — > {7, Z}*(Q') + g{jp$). The phase space element of the second step 
{7, Z}*(Q') -> ^"(ps) + ^+(p 4 ) of the cascade is similar to flT2]): 



(2) 



'47T/X 2 



16tt 2 



ra-e* 







C dy R y R ~ £ (1 - y R )~ e F dcp R . (19) 
jo Jo 



Analogously, we obtained a factorized expression for the hard gluon NLO correction to the 
neutral current process. It has exactly the same structure as Eq. (fl5l) . but the cross section 
in Born approximation &q (zs,e) of the neutral current process has a different form: 



where 
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Integration over y gives for o"H^. d (s, e) the same result as (flEj) . So, in the limit e->0we come 
to the expression almost the same as ([Tj 
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Comparison of this expression where the collinear divergence appears as a pole -, with 

e 

(T5]) where the collinear divergence manifests itself in the form of logarithms permits us to find 
what expression one has to subtract from a^ rd (T5]), namely 
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where m q is the mass of the pair quark and antiquark coming from the both protons. 



2.3 Virtual and soft gluon contribution 

Working in our massive quark treatment we obtained for the sum of virtual and soft gluon 
contributions an expression free from infrared divergences. In the case of charged current 
processes we have 
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Correspondingly, in the case of neutral current processes we have 



a 



NC 
Virt 



NC 
Soft 



a 

71 



S r -NC/ 



+ In 



4u 



-X 



In 



1 



- 2 



1 

2 



- 1 - 



Ti- 
li 



TV 



.. (24) 



(25) 



One can find the collinear divergent expressions to be subtracted from these virtual and 
soft gluon contributions because they are already included into PDF, taking the correspond- 
ing expressions which one has subtract from the hard gluon contributions to the considered 
processes. But one has to take them with opposite sign, to substitute the argument of Born 
cross sections taken z = 1, namely 
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and to integrate over z from to z max . In this way we obtained the expressions to be 
subtracted from virtual and soft gluon contributions. 

For charged current processes collinear divergent subtraction is: 
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And for neutral current processes it is: 
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Integration over z gives the expressions to be found, namely, first - a subtraction for the 
charged current contribution: 
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and second - a subtraction for the neutral current contribution: 
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One can see that subtractions (129]) and (I30I) really subtract the collinear divergent terms in 
the expressions of virtual and soft gluon contributions ff24|) and ff25l) . correspondingly. 

3 Quark— gluon Drell— Yan processes 
3.1 Charged current quark— gluon processes 

In the framework of the Drell- Yan process pp — > W — > Ev£ we have to take into account 
the presence of gluons in the protons. Therefore we consider on the quark-parton level the 
attendant processes with incoming gluon u{p 2 ) + g{p$) — > d(pi) + feips) + i + (p4), see the 
Feynman diagrams on the Figfl], and d(px)-\-g{p§) — > u(p2)+ve(p3)+£ + (p±) (similar diagrams) . 



g{Pb) 





Figure 1: Charged current diagrams with coming gluon. 



The contributions of these processes do not contain infrared divergences but have quark 
mass singularities. 

The three particle phase space element of the process u g — >• d £ + vg can be treated as a 
cascade analogously to (jUJ) but here gluon g is incoming and d - quark is outgoing. The phase 
space element of the first step g{p$) + u{p2) — > W*{Q') + d(pi) of the cascade is the same as 



in Eq. (Tl0|) . Variable z has the same meaning, z 



and s' = —(p3 + Pa) 2 is the invariant 



mass of outgoing leptons, but here s = — 2{p 2 .p^). Variable y has the same form ffTTj) and 9 g 
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is the angle between vectors p\ and p$, but now p*5 is the momentum of the incoming gluon 
and pi is the momentum of the outgoing quark. 

The phase space of the second step W*(Q') — > vzipz) + £ + (Pi) of the cascade, obviously, 
has the same form (|12|) . 

In our massive quarks treatment in the 4-dimensional phase space we obtain the cross 
section of the process u g — > d ug £ + in the form factorized to the Born cross section of the 
corresponding quark- ant iquark process u d — > vg £ + : 
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(31) 



where 



P qg (z)=T f z 2 + (l-z) 



(32) 



is the quark-gluon splitting function and Tf = ~. We neglected the quark masses except 
of in the logarithm, where we have a mass of the outgoing quark, namely we face a mass 

singularity. Integration over z here is up to z max = ( 1 p ] ~ 1. 

V VsJ 

In the MS scheme with massless quarks working in the n-dimensional phase space we 
obtain a factorized expression: 
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We divided here by 1 — e because of the number of the gluon spin projections in the 
n-dimensional phase space is n — 2 = 2(1 — e). After integration over y we have 
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In this expression the collinear divergence appears as a pole -. It has to be compared with 

e 

the analogous expression (}3~Tj) where the collinear divergence appears as a mass singularity. 
In the limit e — > we have 
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Comparing with (13T1) we see that the collinear divergent term which in MS scheme has to be 
subtracted from the cross section of the considered process (because it is already included 
into PDF of gluons) in our treatment has the form 
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Here ma is the mass of the outgoing quark. 

We have the same situation with the process g d 



u v t 



Analogously we obtain that 



from the cross section of this process we have to subtract the corresponding expression: 
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In this process anti-up quark is outgoing and it can move collinearly with the gluon, so the 
logarithm is from the mass m u of the w-quark. 



3.2 Neutral current quark-gluon processes 

In the case of the DY processes pp —>■ [A, Z] —>■ i~ i + we consider on the quark-parton 
level the additional process q{p2) +5 , (ps) — ► q(pi) +^~{p^) +^ + (P4) (corresponding Feynman 
diagrams are similar to the Fig(T]) and process q{p\) +g(p5) — * q(p2) +^~(P3) +C- + {pi) with an 
incoming gluon. The contributions of these processes also do not contain infrared divergences 
but have quark mass singularities. We considered these processes analogously as the charged 
current quark-gluon processes. Difference is only in the Born cross section to which they are 
factorized. 

In our massive quarks treatment in the 4-dimensional phase space the cross section of 
both kind of neutral current quark-gluon processes have equal form: 
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(39) 



where m q is the mass of incoming as well as outgoing quark (or anti-quark). The upper limit 



of integration over z is the same, namely z„ 



1 



1. Here we also neglected the 



quark masses except of under the logarithm. 

In the MS scheme with massless quarks working in the n-dimensional phase space we obtain 
the same expressions fl33l) . fl34|) . fl36l) . as those of charged current quark-gluon processes. 

Analogously comparing gave us that the collinear divergent term which in MS scheme has 
to be subtracted from the cross section of the both considered processes, because is already 
included into PDF of gluons, in our treatment has a form 
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4 Numerical calculations on hadronic level 



4.1 Hadronic level kinematics 

In the c.m.s. of the quark-quark or quark-gluon pair (see Figj2^) we have (neglecting masses 
of quarks): 



Pi = Pi = 1 Pi 1 = 1 P2 



P2 = -Pi, 



s = -2p 1 -p 2 . 



(41) 





a) b) 
Figure 2: a) Quark c.m.s. b) Proton c.m.s. (Born approximation). 



s = -2p lN -p 2N 



In the c.m.s. of protons p,p (see FigfSb) we have 

Pin = P°2N = 1 Pin 1 = 1 &n 1= y) V%n = -&ln> 

Pi=Xip lN , P2 = x 2 p 2N , s = -2pi -p 2 = -2p lN ■p 2N x 1 x 2 = sx!X 2 . (42) 

Equations of transition from quarks c.m.s to protons c.m.s. for any momentum have the form 

Q° = i{Q°n-(3Q n ), Q z = i(Q z n -(3Q° n ), (43) 



where 
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Replacing the scalar products 
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(44) 



(45) 



where 9 q is an angle between the 3-momenta of the quark p\ and neutrino p^ , from quarks 
c.m.s. to protons c.m.s. 

PvP3 = XiPiN'Pz = -XiP° 1N P3° (1 - cos6»jv) , p 2 -p 3 = x 2 p 2N -p 3 = -x 2 p° 2N p z ° (1 + cos6»jv) 

(46) 

and the angle On is between 3-momenta of the proton pi N and neutrino p 3 , after some algebra 
we come to the relation between angles 9 q and On 



cos 0„ 



X 2 — Xi + (xi + x 2 ) cosO 



N 



X X + X 2 + (x 2 - Xi) COS0 N 
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(47) 



The same connection we have for the angle O^ 4 between the quark momentum pi and 
the charged lepton momentum p 4 and for the angle 9]f correspondingly between the proton 
momentum p lN and the charged lepton momentum p 4 . We have only to take into account 
that for the massive charged lepton we have multiply cos^g and correspondingly cos9^ by 

— — , what is not equal to 1 for a massive charged lepton. The same is valid for the neutral 

current case when both final leptons are massive. 

Working with particle momenta components in the protons c.m.s. we can determine the 
parameters needed for the presentation of results on a hadronic level: 

• the transverse momenta of the charged lepton and the missing transverse momenta of 
the neutrino. 

P^n = vW) 2 + (p^n) 2 and p 3 j, = vW) 2 + (P3 V n) 2 - (48) 

• the rapidity of the charged lepton 

1 , Pan + P4r 
7] — - m — n 

2 p 4 % - p 4 z N 

• the transverse invariant mass of the final leptons 

m l = vV^a -cos o, 

where ip 3 ^ is an angle between the final leptons in a plane transversal to z axis 



(49) 
(50) 



4.2 Integration over variables x\ and X2 

We have to integrate over variables x\ and x 2 the whole contribution to the cross section of 
the DY charged current as well as neutral current processes. The leading order contribution 
is 

a LO = Y, [ dx l f(x l ,^ 2 ) f dx 2 f(x 2 ,fi 2 )cTo C (xiX2s), 

V / dx 1 f(x 1 ,ii 2 ) / dx 2 f{x2,n 2 )aQ C {x 1 x 2 s). (51) 
„ Jo Jo 



NC 



In the next to leading order (NLO) we have to add one-loop corrections: hard gluon 
contribution, virtual and soft gluon contribution and corresponding subtractions of quark- 
antiquark processes and also the contribution of quark-gluon processes with their subtractions. 



CX {XiX 2 S) 



a NLO = dx 1 f\x 1 ,{i 2 ) / dx 2 f(x 2 ,H 2 ) 

9192 Jo Jo 
+ ^Hard^l^ S) - <^3 ubt > 2 , X lX2 s) + ^SoftVirt (^2 s) - ^SV Suht M^ X 1 X 2 s) 
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+ V / dxif{xi,^ 2 ) [ dx 5 g(x 5 ,ii 2 
qig Jo Jo 
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(52) 



The subtractions are needed because the MS scheme collinear divergent terms are already 
included into the quark distribution functions f(xi,fi 2 ) and /(x 2 ,/i 2 ) and into the gluon 
distribution function g(x^,fi 2 ). For neutral current processes we have the same formula. 

Let us take into consideration the following expression with quark-antiquark hard and 
soft-virtual subtraction terms: 



[xix 2 s) - ^g ubtr (/i 2 , x x x 2 s) - <x££ Subtr (/i 2 , x x x 2 s) 
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We used the "+" prescription because we have here difference <7q c (z, x±, x 2 , s)— <t^ c (1, Xi, x 2 , s). 
One can note that we can apply subtractions to the PDF instead of the cross section: 



ffa, n ) -> f(xi, v)-7T dzf[—^ 
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2tt 
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The same manipulation can be done for the neutral current processes. In practical applica- 
tions, if the subtraction is applied to PDFs, one has to take care on spurious O (a 2 ) contri- 
bution. The can be done by means of linearization procedure described in Ref. [12J. 

We have a possibility to integrate numerically over three angles and over the independent 
variables z, X\ and x 2 to obtain ^Hard( x i x 2 s), &9^1(xiX 5 s) and o"^(x 2 x 5 s). 



We introduced a new variable W x = X\X 2 when integrated the hard gluon contribution 
of the quark-antiquark process to obtain a^ rd (xix 2 s). If we take x 2 ,Wx for the set of 
independent variables, the Jacobian of the transition is equal to — . So we have the following 
transition of the integrals: 



dxi \ dx 2 f(x 1 )f(x 2 ) 



f l dW x f 1 dx 2 -f(^) f(x 2 ). 
o Jo x 2 \x 2 J 



(55) 

When we integrated the hard gluon contribution of the quark-gluon process to obtain 
&g^(xxX5 s) or frq^ g (x 2 X5 s) we introduced a new variable W y = x\x^z or W y = x 2 x$z. In this 
case we have the following transition of the integrals: 



dxi / dx 5 f(xi)g(x b ) 



dz 



dW„ \ dx 



Zmax ] Wii 

dz — f(^L)g(x 5 ). (56) 

X 5 Z X 5 Z 
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5 Conclusions and Numerical Results 



For the sake of comparison with MCFM [11] for numerical evaluations we used the following set 
of input parameters: 



G f 


= 1.16639 x l(T 5 GeV- 2 , 


a(0) 


= 1/137.03599911, 


oc s {M z ) 


= 0.130, 


a s {M w ) 


= 0.1326, 


M w 


= 80.419 GeV, 




= 2.06 GeV, 


M z 


= 91.188 GeV, 


r z 


= 2.49 GeV, 


M H 


= 115 GeV, 


m t 


= 170.9 GeV, 


m u 


= ma = 66 MeV, 


m c 


= 1.5 GeV, 


m s 


= 150 MeV, 


m b 


= 4.62 GeV, 


\v ud \ 


= |vy = 0.975, 


\v us \ 


= 1^1 = 0.222. 



The CTEQ6L1 [21j set of PDF was used with the factorization scales being equal to M z and 
M w for the NC and CC cases, respectively. The following cuts on the final state kinematics 
were applied: 

P t > 25 GeV, M u > 20 GeV, r) < 1.2, (58) 

where Pt is the transverse momentum of a lepton, Mu is the invariant mass of a charged 
lepton pair (only for NC), and rj is the pseudo-rapidity of a charged lepton. 

In numerical evaluations we used an adaptive Monte Carlo integrator based on the VEGAS 
algorithm [22]. For the partonic sub-process cross sections we used the standard SANC 
FORTRAN modules which can either produced interactively by the system or just downloaded 
from the SANC webpage [4J. These modules are described in Ref. |23j . 

In Fig. [3] we show comparison of the SANC results with the MCFM ones for the transverse 
momentum distribution of fi + in the charged current DY process at LHC. Fig. H] shows the 
corresponding comparison of results for the invariant mass distribution in the NC case. 

Note that the deeps in the first bins of both the distributions are not physical, they appeared 
due to kinematical cuts imposed just at the left borders. Application of the PDF factorization 
in SANC and MCFM are performed in different schemes: the scheme with massive quarks in 
SANC (as described above) and the MS scheme with massless partons in MCFM. We see that 
for the given distributions the difference between these schemes is not numerically important. 

For a realistic application, one has to take into account also QCD showers. It can be 
done with help of the standard packages like PYTHIA [54] and HERWIG [25J. Note that 
the showers will wash out the negatively weighted events, which can be seen in the resonance 
region in Fig. [31 

In this way we presented in detail the evaluation of NLO QCD corrections to Drell- 
Yan like processes. It is important that we performed it in the environment of the SANC 
system, so that now we have a self-consistent simultaneous treatment of QCD and electroweak 
radiative corrections the the DY processes. It is required for the forthcoming experiments 
at the LHC. Simultaneous implementation of the electroweak and NLO QCD corrections to 

2 The LO PDF were used just for the comparison. For practical applications of the described results NLO 
PDF should be chosen. 
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Figure 3: /i + transverse momentum distribution in CC Drell-Yan. 



> 

0) 

o 

2 60 
a 



E 50 
v 




U 



x sand loop 
— mcfm 



i|HililiiilHlHnlri"lll»lT1T"|-|M«1Wltt 




%0 55 60 65 70 75 80 85 90 95 100 

M rr ,GeV 



Figure 4: Lepton pair invariant mass distribution in NC Drell-Yan. 

Drell-Yan processes received with help of the SANC into a Monte Carlo event generator will 
be described elsewhere 
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